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We study the delocalization transition for spin-polarized 
and spin-degenerate non-interacting electrons in the lowest 
Landau level. We perform finite-size scaling calculations of 
the Thouless number, for varying amounts of potential and 
spin-orbit scattering. For spin-polarized electrons, we obtain 
a one-parameter scaling function for the Thouless number 
that fits scaled experimental data for the longitudinal resistiv- 
ity. For spin-degenerate electrons with spin-orbit scattering, 
the Thouless number is peaked away from the band center 
by an amount proportional to the strength of the spin-orbit 
scattering. The universality class of the delocalization transi- 
tion for non-interacting spin-degenerate electrons in the quan- 
tum Hall regime is found to be the same as for spin-polarized 
electrons. We also study the density of states and Thouless 
number for the model of pure spin-orbit scattering studied by 
Hikami, Shirai, and Wegner [Nucl. Phys. B408, 415 (1993)], 
which represents a different universality class. 

PACS numbers: 72.10.Bg, 71.55.Jv, 72.20.My 



1. Introduction 

In the quantum Hall effect changes in the quan- 
tized value of the Hall conductivity Uxy as the magnetic 
field or particle density are varied are due to a series of 
metal-insulator transitions The insulating phase 

occurs when the chemical potential lies in a region of lo- 
calized states. In the insulating regime, Oxy is quantized 
to an integer multiple of /h to a very high degree of 
accuracy over a finite range of field values, and the lon- 
gitudinal conductivity axx vanishes ||] . Near the critical 
points, (Jxx is nonzero and axy no longer remains con- 
stant. The widths of the quasi-metallic regimes become 
smaller as the temperature is reduced |6|. A theoreti- 
cal phase diagram |Q for the integer quantum Hall effect 
(IQHE) has recently been proposed |^ for the transition 
between insulating and quantum Hall regimes. 

Metal-insulator transitions are usually described us- 
ing scaling theories. In the absence of a magnetic field, 
scaling arguments indicate that in two dimensions, all 
states are localized in the thermodynamic limit, leading 
to an insulating state In three dimensions, there can 
exist a range of chemical potential for which the sys- 
tem is a metal, separated by a mobility edge from an 
insulating state. The situation is qualitatively different 



in two dimensions in the presence of a strong perpen- 
dicular magnetic field. Semiclassically, the trajectory of 
a charged particle in a magnetic field can be separated 
into a drift motion along equipotential contours, together 
with small gyrations about the drift motion. In the limit 
of smoothly varying disorder, there is a single energy at 
which equipotential contours percolate through a system. 
This corresponds to an extended state [ p^JTl| . Scaling 
theories of the metal- insulator transition in the quantum 
Hall regime have been developed, consistent with the 
idea that conductivity depends on disorder and chem- 
ical potential [pjp^. For the integer quantum Hall ef- 
fect (IQHE), random potential scattering broadens the 
Landau-level peaks in the density of states. However, 
it is known that the regions of extended states are not 
broadened in energy. There is a discrete set of energies 
{Ec\, corresponding (for weak disorder) to the Landau 
level energies (n-f i)?iijJc, at which the localization length 
^ diverges as ^ ^ \E - E^]-" H. 

The localization exponent v has been measured di- 
rectly in experiments by varying the width of Hall bars 
and determining the scaling behavior of the longitudinal 
resistivity pxx in the metallic regime, at very low temper- 
atures [0. It was found that v w 2.3. The scaling behav- 
ior of the IQHE delocalization transition has also been 
measured by studying the resistivity tensor as a function 
of temperature Q. The reciprocal of the width Ai? of the 
transition region of pxx between Hall plateaux, and the 
maximum slope of the Hall resistivity in the transition 
region, dpxy/dB, both diverge with the same exponent 
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The scaling exponent k is related to the localization ex- 
ponent 1/ by a third exponent p, according to k = p/2i> 
|l3| . The exponent p describes how the inelastic scat- 
tering length lin oc T~p, which determines the effective 
sample size p^ , diverges as a function of temperature. 
For spin-polarized electrons, experiments yield k « 0.42 
Q. When combined with the experimentally obtained 
value V « 2.3 from Ref. |14|, this implies p « 2. The 
value of p has been measured by independent means on 
the same sample in Ref. and p w 2 is also found. 

It should be noted that there are other scenarios for 
explaining the value of k. From the point of view of 
dynamical scaling, k = l/zz^, where z is the dynami- 
cal scaling exponent |17 1^. The value 2 « 1 is consis- 
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tent with Coulomb interactions controlling the dynamics 
0,|9|. Polyakov and Shklovskii have argued that 
hopping transport is responsible for the broadening of 
the Pxx peaks. Sondhi and Kivelson JTst equate temper- 
ature scaling to (imaginary) time or frequency scaling, 
and propose that the result k 1/v is due to the scaling 
behavior of the time scale (rather than the length scale 
lin)- Recently, the microwave frequency dependence of 
Gxx in the IQHE was measured, and the width Ai? of Uxx 
was found to scale with the frequency lo like AB cx w'', 
where 7 = 0.41 ± 0.04 for spin-split peaks ||2^. In this 
work, we explicitly perform length rescaling to obtain 
and do not directly calculate the temperature scaling of 
the conductivity. 

It has also been demonstrated experimentally that 
higher derivatives of pxy scale according to 



an amount 2Ec <C F, where F is the disorder broadening 
of the Landau levels, the localization length has the form 
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where denotes the center of transition region (where 
Pxx has a local maximum), and n = 1,2,3 plf . As- 
suming that Eq. (2) also holds for all higher deriva- 
tives, it follows that in the metallic region, the resistivity 
Apxy = p{Bc) — p{B) near Be is a scaling function of a 
dimensionless length scale, £,/lin] of the ratio of the 
localization length to the effective system size . 

A substantial body of numerical work using finite- 
size scaling [p2|-p5[ shows that, for non-interacting spin- 
polarized electrons, sa 2.3 p^pfl-pS], in agreement with 



experiment |14| and with (non-rigorous) theoretical ar- 
guments that ly = 7/3 |30|. However, for experiments 
with lower mobility samples in smaller magnetic fields, 
the disorder broadening can exceed the Zecman split- 
ting. In this case there can occur a spin-unresolved quan- 
tum Hall transition, with the Hall conductance chang- 
ing by Aaxy = 2e^/h between plateaux. For such ef- 
fectively spin-degenerate electrons, experiments find that 
K = 0.21, which is half the size of the spin-polarized expo- 
nent 1 21 3^. In addition, the exponent of the frequency 
dependence of the width AB of axx is also halved for 
spin-degenerate electrons: AB cx uj'' with 7 — 0.20±0.05 
pO| . There is experimental evidence that p is the same 
for spin-polarized and spin-degenerate electrons in the 
IQHE, which, taken at face value, suggests that i' dou- 
bles in the spin-degenerate case [ p^ . A similar approxi- 
mate doubling of the exponent is seen in network model 
simulations when the localization energy is assumed to 
diverge at only one energy p^ . 

However, Khmelnitskii has argued that for the case 
of two overlapped spin subbands, the extended states 
should split ]3^ , and indeed network model simulations 
of spin-degenerate electrons are fit somewhat better by 
assuming that there is such a splitting . Polyakov and 
Shklovskii |jl^ have conjectured that the for an extended- 
state whose energy is split due to the SO interaction by 
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and that one recovers the usual value for v only very 
close to the spin-split energies ±Ec, at sufficiently low 
temperature. Otherwise, an apparent doubling of is 
obtained. Reference has argued that in the limit of 
very smooth disorder, spin-degenerate electrons diverge 
with the same value of the localization exponent, v « 2.3, 
but at two separate energies. According to these argu- 
ments, the erstwhile doubling of v seen in experiments is 
due to assuming that the localization length diverges at 
a single energy rather than at two nearby energies. 

We shall argue from direct numerical calculations that, 
in agreement with theoretical arguments pof and net- 
work model studies |3^ , the localization length diverges 
at two energies ±Ec, that the localization exponent is 
V « 2.3, and that the value of the universal peak conduc- 
tivity for axx is the same as in the spin-polarized case. 
Our calculation of the Thouless number is done within 
a semirealistic microscopic model, and provides an al- 
ternate confirmation of the energy splitting produced by 
spin-orbit scattering. We have also investigated the pure 
spin-orbit scattering model of Hikami, Shirai, and Weg- 
ner which is believed to be in a different universality 
class. 

2. Computational method 

We have calculated the Thouless number for electrons 
moving in the continuum, restricted to the lowest Lan- 
dau level (LLL). In this calculation, a noninteracting 
two-dimensional electron gas (2DEG) is subject to a per- 
pendicular magnetic field B = VxA, where A is the 
vector potential. The results presented here were ob- 
tained from finite-size calculations in the LLL, using a 
basis of continuum states. For a continuum system with- 
out disorder in the presence of total flux <I>, the density 
of states is a discrete sum of delta functions at energies 
E ^ {n + ^)h.ujc, where n = 0,1,2,... is the Landau 
level index, and lOc = eB /mc is the cyclotron frequency. 
Each Landau level has a degeneracy of — ^/(j>o, where 
(pQ — hc/e is the elementary flux quantum. 

In the Landau gauge A = Bxy, the single-particle 
wave functions are given by 
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where k — 27rm/iV indexes the states within a given Lan- 
dau level, Hn denote Hermite polynomials, and the unit 
length is defined by the magnetic length Ib — ^Jfic/eB. 
The area of the system is measured by the num- 
ber of fiux quanta, N , that pass through the sample, 
— . The Hamiltonian for the system in the 

LLL may be written as 
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where the basis set {|fc,s)} are the n = (LLL) eigen- 
states in Eq. (3), with spin quantum number Sz — fis, 
where s = ±1/2. The projected Hamiltonian in Eq. (4) 
consists of both a random Gaussian scalar potential V 
and a spin-orbit scattering term Hso- The disorder cor- 
relation lengths of V and Hso are taken to be the same, 
and are denoted by C,: 



V{y)V{v') = 
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where U gives the strength of the scalar disorder. The 
limit C ~* corresponds to uncorrelated white-noise dis- 
order. As we discuss in detail further below, we take the 
SO scattering to have the form of a random pseudofield 
W coupled to the electronic spin a: 
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where A is a phenomenological SO coupling constant, and 
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where W gives the strength of the random pseudofield, 
and a,(3 — X, y, z. 

Our numerical studies measure delocalization by cal- 
culating the disorder-averaged Thouless number T{E), 
which has been argued to be proportional to the disorder- 
averaged longitudinal conductivity a^x (E) | |36[ | . The pro- 
portionality of T{E) and (Jxx{E) has been demonstrated, 
in the absence of a magnetic field, in Ref. |^ In prac- 
tice, we express the Hamiltonian as an iV x iV matrix, 
which we diagonalize to obtain iV energy eigenvalues. In 
Ref. 1^^, T{E) is defined as the change AE in the energy 
eigenvalue at energy E that results from changing bound- 
ary conditions {e.g., from periodic to antiperiodic), mul- 
tiplied by the total density of states D{E). The Thouless 
number is small for localized states because such states do 
not extend to the boundaries. It is relatively large for ex- 
tended states because they span the sample to its bound- 
aries and hence are sensitive to changes in the boundary 
conditions. The width Ai? {e.g., at half maximum) of 
T{E) measures the energy range over which states are 
extended. The width is finite for finite-size samples, but 
vanishes in the limit that the sample size becomes infi- 
nite, since the extended states for an infinite sample occur 
at discrete energies, in the presence of a strong magnetic 
field. 

The effect of changing the boundary conditions of wave 
functions is equivalent to the gauge transformation that 
results from adiabatically inserting fiux in a toroidal ge- 
ometry ]3^ . This may be implemented via the perturba- 
tion Hamiltonian 
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where 5 A. is the vector potential for a solenoid carrying 
fiux (p. Insertion of the fiux </> = '/>o/2 is equivalent to 
changing the boundary conditions from periodic to an- 
tiperiodic. In the absence of a magnetic field, the change 
in energy IS.E is second order in 5H ||3^], and /S.E may 
be calculated from second-order perturbation theory in 
5H, which, like the Kubo formula for the conductivity, 
has matrix elements involving the square of the current 
operator. Ando has argued that the the derivation 
given by Licciardello and Thouless |^ relating the Thou- 
less number in the absence of a magnetic field can also 
be applied to the case of a system in a strong magnetic 
field, despite the lack of time reversal symmetry. 

Reference |Q assumed that the mean-square value of 
the matrix elements of the current operator between two 
states is not too sensitive to the difference in energy be- 
tween the states, and that the energy levels are uncorre- 
lated. The latter assumption is not true in the metallic 
regime |^]. Nonetheless, more recent work |37) has ex- 
tended the Thouless formula to the metallic regime, and 
shown that the dissipative (transport) conductance gd is 
proportional to the level curvature gc, defined as 



g,{E) = D{E) 
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In the presence of a magnetic field, it is also believed that 
T{E) and axx{E) are proportional, although we know of 
no proof of this. Differences appear when a magnetic field 
is present; for example, the energy shifts due to changing 
the boundary conditions are linear rather than quadratic 
in (j), due to the lack of time reveral invariance, and are 
random in sign. 

The localization exponent v is obtained by finite-size 
scaling. We calculate the Thouless number for a range of 
samples sizes L oc \/N . In most of our calculations, we 
have followed the approach of Ref. |3^ and calculated the 
energy shifts between periodic and antiperiodic boundary 
conditions, studying 2A£'/0^ rather than d'^E/d4P' , with 
(f) = 0o/2. We expect that T{E) should have the same 
scaling behavior as axx{E). For samples large enough 
to be in the scaling regime, the Thouless number can be 
written in terms of a scaling function of the ratio of the 
localization length ^ to the sample size L, f{^/L). As- 
suming that the localization length diverges at the LLL 
center as ^ oc we may write the Thouless number 

as a scaling function of \E\L^^'^: 

T{E) = f{i/L) ^ mE\L^/^r'^] ^ f{\E\L'^n- (11) 

The area A{L) under the Thouless number curves there- 
fore scales like L^^^^: 

/•oo poo 

A{L)= n{E)dE = L-'/'' f{E)dE=foL-'/'', 
Jo Jq 

(12) 



3 



where T(0) = /(O) and /o are independent of system 
size. T{E) should therefore have a half-width AE cx 
L^^/" that narrows with increasing system size, and a 
peak value r(0), corresponding to the peak conductivity 
o'xx — /2h 1^^, that is independent of the system size 
L. We analyze our results by calculating A{L) for each 
sample size L and we check that we are in the scaling 
regime by ensuring that the Thouless number data all 
scale onto a single curve when the energy E is rescaled to 
E/A{L). The scaling exponent ly can be obtained from 
a log-log plot of A{L) versus L, which, in the scaling 
regime, should approach a straight line with slope —l/i^. 
Studying the scaling behavior of A{L) has the advantage 
that A{L) is less noisy than T{E), since A{L) is obtained 
by integrating T{E). 

3. Spin-polarized electrons 

We first calculate the Thouless number for spin- 
polarized electrons in the LLL, for white-noise disor- 
der (C = 0). The Thouless number plots for N — 
20, 80, 300, 1000 are shown in Fig. la. The half-width 
AE of the Thouless number shrinks with increasing sys- 
tem size, while the peak value remains constant at about 



0.32 « I/tt. Ando's formula, Gx 



{e^/An)T{E), 



lating the Thouless number to the conductivity for the 
case of a 2DEG in a strong magnetic field |2^, would 
imply that T(0) w I/tt corresponds to a peak conduc- 
tivity (T^^ « e^/2/i, consistent with the known value for 
this model |^^. However, it should be noted that the 
value of the peak Thouless number depends on the de- 
tails of the method used to calculate the energy shift, and 
varies slightly with the correlation length ^ of the disor- 
der. Nevertheless, we expect the peak value of T{E) to be 
independent of system size, and the shape of the scaling 
function and the value of the localization exponent v to 
be independent of the precise way the Thouless number 
is defined. 

We have checked that T{E) is in the scaling regime by 
rescaling the Thouless number plots for different system 
sizes onto a single curve, shown in Fig. lb, for seven sys- 
tem sizes, from = 20 to 1000 flux quanta. The data 
is collapsed onto a single curve by rescaling the energy 
according to ^ E/A{L). A log-log plot of A{L) versus 
L is given in Fig. Ic. A least-squares fit to a straight 
line gives a slope that corresponds to v = 2.36 ± 0.02 for 
the spin-polarized LLL with white-noise disorder. How- 
ever, closer analysis shows that the data deviate from a 
straight line, presumably due to finite-size corrections to 
scaling and so we believe that the error estimate 
for our numerically obtained value of v should signifi- 
cantly larger than 1%. In any case, the value we ob- 
tain of V w 2.4 is consistent with previous numerical 
work that used different methods and sample geometries 
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We have compared our numerical data for the Thou- 
less number to experimental longitudinal resistivity data, 
taken from data used in Ref. [Q. Figure 2a shows the 



experimental longitudinal resistivity data versus applied 
magnetic field, for an InGaAs/InP heterostructure at 
four temperatures, T = 40, 110, 305, 640mK. The sym- 
bols in Fig. 2b are the data of Fig. 2a, centered, nor- 
malized by the peak resistivity, and rescaled, for four 
different temperatures |^^. The width of Pxx{B) nar- 
rows as the temperature becomes smaller. It does so as a 
power law of the temperature, with associated exponent 
K K, 0.42 ± 0.04. Under rescaling of the magnetic field 
with the temperature according to AB AB ■ T~'^, all 
the longitudinal resistivity curves can be collapsed onto 
a single curve. 

Figure 2b shows that the scaled Thouless number 
data fits the scaled resistivity data. The fit between 
the temperature-rescaled pxx data and the size-rescaled 
T{E) curves is further evidence that the temperature 
and length are related by an inelastic scattering expo- 
nent value of p — 2, which has also been found in other 
experiments |l6| . Figure 2b is the universal crossover 
function describing the transition from metal to insula- 
tor. Reference also used the network model to calcu- 
late a crossover function as a function of the localization 
length ^, and obtained good agreement with the exper- 
iments of Ref. Q], after correcting for the contribution 
due to edge currents. Experiments of Hughes et al. ]4^ ] 
find single-parameter scaling of axxiB, T), with k « 0.42 
and axx well fit by the model of Ref. |^ . The experimen- 
tal data in Fig. 2 were obtained from samples of width 
600 microns; we expect edge currents in these samples to 
be less important than for the samples used in Ref. [Q , 
which were only 7 microns wide. If Coulomb interac- 
tions dominate, one estimate for the sample size at which 
edge effects become important is when the sample width 
is smaller than the electron scattering length Zy, where 
e'^/eh ^ ksT. For T=20mK, this gives h ^ SO/xm. 

4. Spin-unresolved quantum Hall transition 

Reference [jsil studied the scaling of the maximum 
slope of the Hall resistivity as a function of temperature, 
for spin-unresolved quantum Hall transitions. Recall that 
for the spin-resolved Landau levels the scaling exponent 
IS K — 0.42 when the spin-splitting is not resolved, it 
is found that k — 0.21 @,^ltl. This would correspond to 
doubling the exponent v sa 2.3 to v w 4.6, assuming that 
the inelastic scattering length exponent p is unchanged. 
Other experiments [p^ give evidence that p = 2 for both 
the spin-resolved and spin-unresolved cases. 

In practice, electrons can be spin-flip scattered, even in 
the absence of magnetic impurities, due to the spin-orbit 
(SO) interaction. The SO interaction is a relativistic ef- 
fect, and is therefore small in comparison to the random 
potential scattering. Nonetheless such scattering could 
be important when the spin levels are degenerate, and 
is stronger for materials with scatterers of high atomic 
number, such as the InGaAs samples used in Ref. [ pl| . 

SO scattering is due to the fact that an electron mov- 
ing with velocity v in the presence of an electric field 
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E = i experiences in its rest frame an effective mag- 
netic field B = ivxE which couples to the electron spin 
according to (e?i/mc)B • (§)<t. The orbital motion of an 
electron therefore gives rise to an effective magnetic field 
that couples to the spin. The SO interaction is written 
as 



^so = 



cr • VVXU, 



(13) 



where 11 = —iK^ + (e/c)A is the mechanical momentum 
of the electron, whose spin is altered by the effect of the 
Pauli matrices a in Eq. (13). 

The effect of Vso on the Bloch electrons of bulk III-V 
compounds such as GaAs has been investigated in Refs. 

In the vicinity of the F point, the effective con- 
duction band Hamiltonian becomes 
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where m* is the conduction band effective mass, Kx = 
k^(k1 — ky), and Ky and Kz are obtained from the ex- 
pression for Kx by cyclic permutation. In a material that 
lacks inversion symmetry, the SO interaction [Eq. (13)] 
leads to a nonzero value of the constant /3, and hence to 
a splitting of the conduction band states which is pro- 
portional to the cube of the wave vector. Application of 
Eq. (14) to inversion layers, developed in Refs. p8| , p9[ , 
leads to the following effective SO Hamiltonian: 



2m" 



+ aia'^py - (T>x) + Picr'^Px ~ (J^Py) (15) 



The third term is identified as P{cr ■ k), the expecta- 
tion value of the effective conduction-band SO interac- 
tion taken in the appropriate quantum-well bound state 
(the interface is assumed normal to z). Khaetskii | |4^ ] 
finds /3 ~ 1.1 X 10^ cm/s for an inversion layer density of 
n — 2 X 10^^ cm^^. The second term arises due to the 
asymmetry of the confining potential, which generates an 
electric field in the inversion layer. Only rough estimates 
of the constant a exist [^,|^, with a assumed to be in 
the range 0.1/3 to /3. In any case, for typical momenta 
of size h/ls, where Ib ~ lOOA, the SO energy scale is 
roughly 10"^ eV, or 0.1 - 1.0 K. 

The SO Hamiltonian in Eq. (15) produces random 
scattering of the spin in the presence of ordinary random 
scalar potential scattering. This is because scalar scat- 
tering changes the momentum of the scattered electron, 
which by Eq. (15) is coupled to its spin. We estimate 
direct SO scattering from the slowly varying impurity 
potentials in high mobility GaAs heterojunctions to be 
extremely small and unlikely to have any significant ef- 
fect |51]. However, in the InGaAs samples studied in 
Ref. (which have a mobihty of 35,000 cmVVs), alloy 
scattering from short-range potential disorder produces 



larger momentum transfers and hence larger spin-flip am- 
plitudes. In addition, the relatively large atomic number 
of In may further enhance SO scattering. Nevertheless, 
the SO energy scale is much smaller than that for scalar 
potential disorder. 

In the presence of an external magnetic field B = 
VxA — Bz, one must add a Zeeman term Hz = 
— ^gfisBa^ , as weU as make the substitution p ^ 11 = 
~ihW + (e/c)A. We have included electron spin in our 
Thouless number calculations, in the limit of zero Zee- 
man splitting ((7 = 0), so that the spins are not resolved. 
Since 11 is not a constant of motion for an electron mov- 
ing along an equipotential contour of the random po- 
tential V , we take, for simplicity, the SO coupling to 
have the form of Eq. (7), a Zeeman coupling to a random 
pseudofield W. In the absence of spin-orbit scattering, 
the exponent v remains the same as for the spin-resolved 
case, although T{E) doubles, reflecting the doubling of 
the density of states due to the spin degeneracy. 

The use of Eq. (7) as the SO interaction rather than 
the more cumbersome Eq. (15) can be motivated by con- 
sidering the Hamiltonian H = Hq + V, where V is a 
random scalar potential, and 
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This is just Eq. (15) in the presence of a magnetic field, 
and in the absence of the term proportional to a. (For 
convenience, we have performed a tt rotation about the 
X-axis in spin space to change —ay ay in Eq. (15).) 
Hq may be written as 

Ho = hijJc{a)a + 1/2) + -^(a^cr^ + aa+), (17) 

where the lowering operator is defined as 
Ib 



V2n 



(H, - iUy), 



(18) 



the raising operator is the hermitian conjugate of a, and 
a± = i iay . Hq can be diagonalized exactly, and has 
eigenstates of the form 



\k+)=u+\2k,^)+v+\2k + l,i) 
|fc-) =7/„|2fc,i)+t;_|2fc-l,t). 



(19) 



where |n, s) denotes a state with Landau level index n 
and spin s (s =t, i). We have omitted the indices for the 
iV-fold degenerate states within each Landau level. 

We now restrict our consideration to the energy eigen- 
states of Hq with the two lowest energies. The other 
eigenstates are higher in energy by multiples of hujc ■ The 
lowest energy states are 



l+>-|0,i), 
|-) = (l-7V2)|0,T)-7|l,i), E_ 



= 



(20) 
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where 7 = (3 / \/2lBi^c ^ 10~^, and we measure the ener- 
gies E± with respect to huJcf^. 

In the subspace of these states, it fohows from Eq. (19) 
that the effective scattering Hamiltonian can be rewrit- 
ten, up to a constant term, as 

H,ff ^V + \W -(T, (21) 

where V is just the scalar potential V projected onto the 
LLL, and, to order 7, 

iW^,Wy) = ~27(Rc(+|V|-),Ini(+|V|-)). (22) 

Wz is of order 7^, but is not expected to change the 
universality class of the delocalization transition in the 
presence of the scalar potential V which, like a Wz^z 
term, produces scattering without flipping the spin. 

If SO scattering were a relevant perturbation (one that 
changes the universality class of the transition), its effects 
would become large at sufficiently long length scales and 
low enough temperatures. We check the relevance of SO 
scattering by artificially increasing the strength of the SO 
scattering, and then calculating the exponent v and the 
peak value of the Thouless number r(0). When short- 
range SO scattering is included, we find that the half- 
width A_B of the Thouless number curve becomes larger 
and has a scaling exponent of about v = 4.4±0.2, assum- 
ing that the localization length diverges only at the band 
center, as for the spin-resolved case. This is illustrated in 
Fig. 3a by the Thouless number curves for system sizes 
of = 20, 80, and 500, which are peaked at = 0. Fig- 
ure 3b shows a log-log plot of AiV) versus L, as in Fig. Ic. 
The apparent value of v is twice that of the spin-resolved 
value, in agreement with the experimental results of Ref. 
pl| . The peak value of the Thouless number, T(0), in- 
creases from 0.32 to roughly 0.4. This is similar to the 
behavior seen in the experiments of Ref. jsi] and in the 
numerical simulations of Ref. ||3^ , when ^ is assumed to 
diverge at a single energy. This would seem to suggest 
that SO scattering might belong to a different universal- 
ity class than ordinary scalar potential scattering. 

However, when the SO scattering potential is smoothed 
(C — 2), we obtain very different results. This is illus- 
trated in Fig. 4a, where we have plotted the Thouless 
number for system sizes of iV = 40, 160, and 500. The 
Thouless number curves acquire two symmetrically dis- 
placed peaks at zbi^c (we show only one, for > 0), 
each with a peak value of approximately I/tt, just as 
for the spin-resolved case in the absence of SO scatter- 
ing. If we assume that the Thouless number is the sum 
of symmetrically displaced scaling functions, then a log- 
log plot of A{L) versus L, shown in Figure 4b, yields 
v Ri 2.5, much closer to the value for the spin-resolved 
exponent, ly — 2.3. The similarity of the peak height 
T(0) and the exponent v to the case of spinless electrons 
are evidence that the universality class of the delocaliza- 
tion transition is the same for the spin-unresolved and 
spinless non-interacting electrons. 



We have also varied the strength A of the SO scat- 
tering to measure its effect on the energy Ec at which 
T{E) peaks. This is shown in Fig. 5a, where we have 
plotted T{E) versus E fov N = 80 and C = 2, with 
A = 0,2,4,6. We find that for this simple form of SO 
scattering, Ec oc (A|W|), just as for Zeeman splitting. 
Figure 5a can be understood in terms of an effective 
Zeeman coupling generated by the smooth SO scattering 
pl[ . Denoting the correlation length of the SO scatter- 
ing by and the average strength of the random field 
by (|W|), it can be shown that when ( ^ Ib, so that 
W is slowly varying, the SO field W acts like a local 
constant magnetic field coupling to the spins and gener- 
ates an effective Zeeman splitting that depends on the 
magnitude (but not the direction) of W. The slow vari- 
ation in the direction of W yields a Berry's phase that 
acts like a weak random magnetic field, but this does not 
change the universality class in the presence of a strong 
external magnetic field j2^j5^. The net effect, there- 
fore, is to produce a splitting between ±Ec of magnitude 
AE = 2(A|W|). This behavior is shown in Figs. 5b and 
6b. 

It has been suggested that in some IQHE models, the 
universality class of the metal-insulator transition for 
white-noise disorder = 0) is different than for smooth 
disorder ||2|]. In order to investigate whether this is the 
case here, we have artificially increased the value of A 
as in Fig. 5a, except that we use white-noise disorder. 
The results for N = 20, with A = 0, 2, 4, 6 are shown in 
Fig. 6a. It can be seen that, as in the case of smooth 
disorder (Fig. 5a), the peak of the Thouless number is 
displaced away from zero {Ec 7^ 0). A log-log plot of 
A{L) versus L, however, yields v = 4.4 ± 0.2. We believe 
this is due to the slow crossover to the scaling regime 
when C/Ib < 1, as found by Huckestein (26), but cannot 
demonstrate this within the range of available sample 
sizes. 

Huckestein has studied numerically finite-size correc- 
tions to scaling in the center of Landau levels |^^. Al- 
though the value of i/ « 2.3 is obtained for the lowest 
{n = 0) LL, independent of the correlation length C, of 
the disorder potential, this is not the case for the sec- 
ond (n = 1) LL. For n = 1, w 2.3 is observed only 
when C > Ib- For C < no universal scaling behav- 
ior was observed. Following Chalker and Eastmond 
Huckestein showed that deviations from finite-size scal- 
ing themselves scaled with size, so that, in effect, Eq. (10) 
is modified to 

T{E)^ !{\E\L''-,krrL-y-), (23) 

where Urr is a length scale that is a function of C/^s- The 
exponent yirr « 0.38 ± 0.04 obtained by Huckestein p^ ] 
agrees with the value obtained in Ref. for network 
model simulations. Of great practical importance is the 
fact that the size of the length scale Urr is 10"* times 
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longer for C = than for ( = 0.81b- This makes it 
difficult to observe scaling behavior for n = 1 when C = 
0. It is possible that the apparent dependence of v on 
C that we find in the spin-degenerate case may have a 
similar origin to the difficulties encountered in reaching 
the scaling regime when ri = 1. 
5. HSW model 

Hikami, Shirai, and Wegner (HSW) have consid- 
ered a special model for Anderson localization of spin- 
degenerate non- interacting electrons . They study the 
case when only scattering between one spin state and a 
different spin state is allowed. This corresponds to taking 
the pseudofield W to lie in the xy plane, and completely 
supressing ordinary scalar potential scattering. Although 
it excludes important scattering processes present in real 
samples, it is the simplest model of spin-orbit scattering 
in the quantum Hall regime, and has the great advantage 
of being analytically tractable for the case of a Gaussian 
white-noise random distribution of scatterers |]35[| . 

The model of HSW exhibits novel behavior near the 
band center {E — 0). The conductivity at i? = is given 
exactly by b^/tt/i, and the DOS diverges at the band cen- 
ter [HIQ . HSW used a 1/5 expansion, where S is essen- 
tially the electron spin (physically, S" = 1/2), to calculate 
the DOS to order 1/5*^. They found al/S'^ enhancement 
of the DOS at £; = of the form (1/5*2) ln2(|£;/T/RMs|)- 
Using an effective field theory describing the IQHE and 
projecting on to the LLL, Gade finds that the DOS 
diverges instead as a power law, D{E) cx near 
E = {]. We also note that Ludwig et al. ||5^ have studied 
a class of IQHE models which respect symmetries anal- 
ogous to those of the HSW model, and that they also 
obtain a DOS that diverges as a power law near E = Q. 

We have calculated the DOS and Thouless number for 
the HSW model, and we find that both quantities van- 
ish at £* = 0, for finite-size samples. Figure 7a shows 
the DOS of the HSW model as a function of energy 
for the small system size of A^=10 flux quanta. Sys- 
tem sizes of A^=20,40,80,160 also show the same behav- 
ior; namely, a DOS with two strong peaks near E = Eq 
due to the energy level spacing in a finite sample, where 
Eq oc Vrms/A^i corresponding to the average level spac- 
ing between energy eigenvalues. In addition, the DOS 
vanishes linearly a,s E ^ 0, ior E < Eq. We find that 
the peak DOS at £* = £'o increases with N. The HSW 
Hamiltonian may be written in the form 



HhSW = WxCTx 



Wyffy, 



(24) 



where Wa denotes the components of the pseudofield W, 
projected onto the LLL, and the CTq are Pauli matrices. 
Note that crz{HHsw)o'z — ~Hhsw, so that for every 
eigenstate |'0) of energy E, azltp) is also an eigenvalue, 
with energy —E |^^. The energy spectrum of Hhsw 
therefore possesses a reflection symmetry about E — 0. 
The reduction of the DOS near £* = is due to the com- 
bined effects of the reflection symmetry of the eigenvalue 



spectrum, plus the usual level repulsion of eigenvalues. 
It is therefore a finite-size effect. We note that when the 
disorder is smoothed (C > 0), the DOS becomes more 
strongly peaked near E = 0, and that the downturn 
moves closer to £^ = {Eq becomes smaller). This is 
illustrated in Fig. 7b, where we have plotted the DOS for 
a samples of size = 80, for disorder correlation lengths 
C=0,l,2 (in units of the magnetic length). The case of 
smooth disorder was studied in Ref. |Q . 

In order to observe a divergence in the DOS near 
E — 0, one must study energies small enough that 
the HSW correction can be observed, but large enough 
(E > Vrms/N) that the vanishing of the DOS due to 
eigenvalue repulsion does not interfere. Given our nu- 
merical limitations on the maximum size {Nmax ~ 500) 
of our samples, we were unable to isolate any divergent 
contribution to the DOS at small E. To better under- 
stand the behavior of the HSW model near E = 0, we 
studied a random matrix problem with the same sym- 
metry as the HSW Hamiltonian, namely the eigenvalue 
spectrum for matrices with the same form as the HSW 
Hamiltonian, 



H 



f W 

yw^ 



(25) 



where is a Gaussian random complex matrix [ p7[ . 
When the elements Wij all have the same variance, so 
that 

P(Ty,j)«e-l^-l'/2'^', (26) 

one obtains the exponential ensemble of random matri- 
ces, studied by Bronk |^^. The joint probability distri- 
bution function for the eigenvalues is fs^] 



p(i?i, ...,en) (X n i^^.i n - ^kfe 

i j<k 



2 _ k2|2 -X;.-E?/2<T^ 



(27) 



Using a mathematical identity for sums of the product 
of two Laguerre polynomials , one can show that the 
DOS for the exponential ensemble is given by 



DiEi)^ J dE2...J dENP{Ei,...,EN) 
(X \E,\ \l'^},\{x)L%1,{x) L%\{x)Lf{x) 



(28) 



where x — E\l2a'^ and the lI"-* denote associated La- 
guerre polynomials. The DOS is proportional to \E\\ 
near E\ = 0, and vanishes like e^^^l"^" for large 
which is the behavior seen in Fig. 7a for the HSW 
DOS. However, Bronk's calculation of the asymp- 
totic (A^ — > cxd) random-matrix DOS may be invoked to 
show that the usual semicircle law for the DOS holds, so 
that the random matrix model gives no divergent contri- 
butions to the DOS. Not surprisingly, the HSW model 
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and the exponential ensemble have very different behav- 
ior even near E = 0, despite the similar symmetries of 
the Hamiltonians. Random matrix theory is expected to 
apply in metallic phases but not at critical points [59| ] . 

Figure 7c shows the Thouless number for the 
HSW model as a function of energy for system sizes 
A^=10,40,160. It can be seen from Fig. 7c that there 
are two delocalization transitions, one at Ec ~ 0.7Vrms, 
and the other at -E = 0. The transition at Ec ^ OJ^rms 
is reminiscent of the transition at E ^ A|W| studied 
in Sec. 4, for strong SO scattering. As the system size 
increases, it is clearly seen that the peak in T{E) at 
E = Ec 0.7Vrms becomes narrower and more pro- 
nounced. We also find that the peak value, T{Ec), for the 
HSW model is close (only 10% smaller) to that for spin- 
less electrons, T{0). The transition at i5 = is more dif- 
ficult to characterize, since the width of the peak in T{E) 
near E = seems to scale like the energy eigenvalue spac- 
ing, i.e., AE oc 1/N oc L"^, rather than AE cx L^^/'^ 
(near Ec). The delocalization transition near E — Q also 
differs from the one at Ec in that the peak Thouless num- 
ber is seen to (slowly) decrease with increasing system 
size. These differences in the behavior of T{E) are evi- 
dence that the universality class of the two delocalization 
transitions in the HSW model are different. 

6. Conclusions 

Using a continuum LLL model, we find, in agreement 
with previous numerical work, that for spin-polarized 
electrons in the LLL, the localization exponent has the 
value ly « 2.4, and that the peak value of the Thouless 
number is constant, independent of the sample size. We 
find that the scaled Thouless-number data fits the scaled 
experimental resistivity data of Ref. which is the 

universal crossover function for the quantum Hall metal- 
insulator transition. 

For short-range SO scattering, we find, as in the ex- 
periments of Ref. a seemingly doubled exponent, 
v « 4.4, when we assume that the Thouless number is 
peaked at a single energy. This value of v reverts to close 
to the usual non-SO value of v fa 2.3, and the peak value 
of the Thouless number approximates the spin-polarized 
value, when the SO scattering potential is smoothed and 
made sufficiently large. We conclude that the universality 
class of the metal-insulator transition for non-interacting 
spin-degenerate electrons in the quantum Hall regime is 
the same as for spin-polarized electrons. 

It is possible that the experiments leading to an appar- 
ently different value of the temperature exponent k may 
not be in the scaling regime, especially for nearly white- 
noise scattering, which is the situation for the In-doped 
GaAs samples used in Ref. . The experiments of Ref. 
pit required temperatures T > 50mK; otherwise, the 
spin-splitting between the Landau levels was resolved, 
and the putative spin-degeneracy was lost. It could be 
that the effective sample size Im oc T~p was too short to 
reach the scaling regime, even at the low temperatures 



used in the experiments. 

It is not known whether Coulomb interactions change 
the value of k for spin-degenerate electrons, e.g., by 
changing the inelastic scattering length exponent p (or 
the dynamical exponent z) from its value for non- 
interacting electrons iQ. It is believed that for dirty 
boson systems, the dynamical exponent depends on the 
form of the interparticle potential [ |6l| . Certainly a peak 
value of the longitudinal conductivity less than the uni- 
versal value tr^j. = e^ /2h for non-interacting electrons is 
seen in experiments, and perhaps this is due to Coulomb 
interactions. The importance of Coulomb interactions 
even for the spin-polarized case is underscored by the 
large exchange enhancement of the elecron g factor rel- 
ative to its bare value, which has been studied in nu- 
merous experiments |62| and calculations [|63| . Recent 
work has developed a new picture of the nature of spin 
excitations in the quantum Hall regime, in the presence 
of Coulomb interactions |Q. In principle, the 1/r na- 
ture of the Coulomb interaction in a 2DEG could be 
purposely modified by specially gating a heterojunction 
with a metallic overlayer, so as to screen the Coulomb in- 
teraction and induce dipole-dipole interparticle in- 
teractions. This has been carried out for the insulating 
regime, where hopping dominates the transport ^5 66 1. 

We also studied the model of Hikami, Shirai, and Weg- 
ner and calculated the DOS and Thouless number 
for a range of system sizes. We found that the DOS 
and T{E) both vanish linearly at i? = for finite-size 
samples, contrary to the naive expectation that the DOS 
should increase and T{E) become constant at £^ = 0. 
The different scaling behavior of the Thouless number at 
the two energies are evidence that the two delocalization 
transitions are in different universality classes. 
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FIG. 1. (a) Thouless number data for spin-polarized elec- 
trons in the lowest Landau level, in the presence of white-noise 
disorder (C = 0), for sample sizes of 7V=20, 80, 300, 1000 flux 
quanta. The width of the Thouless number decreases with 
increasing A''. 

(b) The scaled Thouless number data for spin-polarized elec- 
trons in the lowest Landau level, for ^^=20, 40, 80, 160, 300, 
600, 1000 flux quanta, falls onto a single curve. 

(c) Log-log plot of the integrated area under the Thouless 
number curves, versus the system size {L oc %/iV), for A''=20, 
40, 80, 160, 300, 600, 1000 flux quanta. The slope of the 
resulting line is and yields v « 2.4. 



FIG. 6. (a) Thouless number data for strong spin-orbit 
scattering (A{|W|) ~ Vrms) and white-noise disorder (^ = 0), 
for spin-orbit coupling values of A=0, 2, 4, 6. 
(b) Ec versus A, for A=0, 1, 2, 3, 4, 5, 6. 

FIG. 7. (a) Density of states D{E) for pure spin-orbit 
scattering model of Hikame, Shirai, and Wegner (HSW) 
[ ^ , for system size A''=10 flux quanta. For finite N, 
D{E = 0) = 0. D{E) exhibits a strong peak at an energy 
-Bo oc 1/A'', close to E = 0, with a peak height that grows 
slowly with A'^. 

(b) Density of states for HSW model, for fixed A=80, and 
varying disorder correlation lengths Cj" = 0,1,2. The peak 
DOS near E = increases with (. 

(c) Thouless number T{E) for HSW model, for system sizes 
N — 10, 40, 160, showing two delocalization transitions, at 
E ~ 0.7Vrms and E ^ 0. For finite N, T{E = 0) = 0. The 
Thouless number decreases with increasing A, except near 
the fixed points E « ±0.7Vrms. 



FIG. 2. (a) Longitudinal resistivity {pxx) data versus ap- 
plied magnetic field for four temperatures, T=40, 110, 305, 
640 mK, taken from Ref. 

(b) Scaled, normalized experimental data for p^x for 
T = 40, 110, 305, 640 mK, shown by symbols. The solid lines 
are the scaled Thouless number data for spin-polarized elec- 
trons in the lowest Landau level, for A=20, 40, 80, 160, 300, 
600, 1000 flux quanta. This curve represents the universal 
crossover function from metal to insulator for non-interacting 
electrons in the quantum Hall regime. 



FIG. 3. (a) Thouless number data with spin-orbit scat- 
tering, white-noise disorder ((" = 0), and sample sizes of 
A'=20,80,500 flux quanta. The width of the Thouless number 
decreases with increasing TV. 

(b) Log-log plot of the integrated area under the Thouless 
number curves, versus the system size. The slope of the re- 
sulting line yields ~ 4.4. 



FIG. 4. (a) Thouless number data for smooth spin-orbit 
scattering, for smooth [C, = 2) disorder, and sample sizes of 
N=AQ, 160, 500 flux quanta. Note that the data are peaked 
away from zero energy, at the same energy Ec. The width of 
the Thouless number decreases with increasing A', 
(b) Log-log plot of the integrated area under the Thouless 
number curves, versus the system size. The slope of the re- 
sulting line yields u ~ 2.5. 



FIG. 5. (a) Thouless number data for smooth spin-orbit 
scattering ((" = 2), for spin-orbit coupling values of A=0, 2, 4, 
6. The Thouless number is peaked at Ec oc (A|W|). 
(b) Ec versus A, for A=0, 1, 2, 3, 4, 5, 6. 
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